Conductance Correlations Near Integer Quantum Hall Transitions 
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In a disordered mesoscopic system, the typical spacing between the peaks and the valleys of the 
conductance as a function of Fermi energy Ef is called the conductance energy correlation range 
Ec. Under the ergodic hypothesis, the latter is determined by the half- width of the ensemble 
averaged conductance correlation function: _F = {5g{EF)Sg{EF + In ordinary diffusive 

metals, Ec ~ D/L^, where D is the diffusion constant and L is the linear dimension of the phase- 
coherent sample. However, near a quantum phase transition driven by the location of the Formi 
energy Ef, the above picture breaks down. As an example of the latter, we study, for the first time, 
the conductance correlations near the integer quantum Hall transitions of which Ef is a critical 
coupling constant. We point out that the behavior of F is determined by the interplay between the 
static and the dynamic properties of the critical phenomena. 



PACS numbers: 73.50.Jt, 05.30.-d, 74.20.-z 

The quantum interference effects in disordered phase- 
coherent systems belong to the mesoscopic physics [0,^. 
A phase-coherent sample is one in which the phase- 
coherence length is larger than the sample size L. 
Thus mesoscopic physics naturally appears in small sys- 
tems of mesoscopic dimensions, usually in nanostruc- 
tures. Mesoscopic physics is also important in systems 
large enough to exhibit macroscopic quantum phase tran- 
sitions. The reason is that L^, being the cutoff for the 
critical fluctuations, diverges as the temperature T ap- 
proaches zero at such transitions. An important ex- 
ample of the latter is a two dimensional electron gas 
(2DEG) close to the transition between two quantized 
Hall plateaus. The mesoscopic fluctuations of the con- 
ductance in this case has been studied recently both ex- 
perimentally and theoretically j^-^. 

The essential physics in the mesoscopic regime is the 
lack of self-averaging in the transport properties. Sam- 
ple specific, reproducible fluctuations in the conductance 
become observable at low temperatures. From a theo- 
retical point of view, the basic statistical properties of 
conductance fluctuations are determined by the conduc- 
tance correlation function H, 



F{Ef,AE) =< SgiEp + AE)6g{EF) >, 



(1) 



where Sg is the deviation from the impurity averaged 
conductance, i.e. Sg{Ep) = g{Ep)— < g{Ep) >. Note 
that the conductances in F are to be evaluated at two 
different Fermi energies separated by the amount AE. In 
general, F can have another argument AB, representing 
the magnetic field correlation of the conductance, which 
we shall not consider here. 

For AE = 0, Eq. (|^) gives the variance of the con- 
ductance, Var((7) = F{0) =< Sg^ >. In the rest of the 
paper, we measure the conductance in units of e^//i and 
the variance in units of (e^//i)^. Under the ergodic hy- 
pothesis, the sample to sample fluctuations are analogous 
to the fluctuations of the conductance as a function of 
the Fermi energy. In this case, the typical spacing of 



the peaks and valleys in the conductance as a function 
of energy in a specific sample, usually called the energy 
correlation range E^ is determined by the half- width of 
the ensemble averaged conductance correlation function 
F{AE) in Eq. {§, i.e. E, = AE^/^. 

In ordinary mesoscopic disordered metals in the dif- 
fusive regime, the disorder-averaged conductance, (g), 
can vary by orders of magnitude, but the variance of 
the conductance assumes a universal value of order one 

. Moreover, the energy correlation range of the conduc- 
tance is given by Ec ~ hn^D/L^, where D is the diffusion 
constant and L is the linear dimension of the system. It 
is important to emphasize that in the diffusive regime, 
the conductance correlation in Eq. (|l|) is only a function 
of the energy difference AE and is independent of Ep 
Ec in this case corresponds to the inverse diffusion time 
across the sample in the current direction which is often 
referred to as the Thouless energy. This is the character- 
istic Fermi energy difference beyond which the paths of 
two injected electrons seize to be phase coherent, giving 
rise to significant difference in the conductances. 

However, as we shall show in this paper, in the criti- 
cal regime of a quantum phase transition (QPT) that is 
driven by the location of the Fermi energy instead of cor- 
relation strengths, the above picture breaks down. The 
primary reason is that, in this case, the Fermi energy is a 
critical coupling constant that controls the proximity to a 
quantum critical point. As a result, the conductance cor- 
relation function in Eq. ([l]) is determined by the critical 
properties associated with the QPT. We shall focus on 
the QHE in which the transitions between the quantized 
Hall plateaus as a function of the magnetic field is driven 
by the location of the Fermi energy of the disordered 2D 
electron system. 

It is well known that the quantum Hall transition 
(QHT) is a continuous zero temperature phase transi- 
tion at a single extended state energy E* between two 
adjacent Hall plateaus [|o|. The critical singularity at 
the QHT is described by a single divergent length, the 



1 



localization length, ^{Ep) cx \Ef — E*\~^ , as Ep ap- 
proaches E* . Here v is the localization length exponent. 
We shall focus on the critical regime, finite systems of lin- 
ear dimension L, and the zero temperature limit. In this 
case, the critical fluctuations are cutoff by L. The width 
of the critical regime shrinks with increasing L accord- 
ing to L"^!'^ . The transition thus acquires a finite width 
W oc L~^l^ . The conductance in the transition regime 
is dominated by phase coherent transport and thus ex- 
hibits mesoscopic phenomena. What is different from or- 
dinary diffusive metals is the proximity to the quantum 
critical point. If the typical spacings between the peaks 
and valleys in g{Ep), i.e. the energy correlation range 
Ec a large number of oscillations would appear 

within the critical region W so long as a > l/i^ ~ 0.42, 
which will be shown to be the case below. In this regime, 
the ergodic hypothesis, which is expected to fail in the 
plateau phases, remains valid. Moreover, the conduc- 
tance correlations are determined by the critical proper- 
ties associated with the QHT. 

We now proceed to write down the scaling form for the 
correlation function defined in Eq. (nh near the QHT, 



F{Ep,AE) = T 



L 



L 



L 



^Ep + AE)' aEp)' L^{AE) 



(2) 



Here is the length scale introduced by a finite fre- 
quency. The origin of the latter is the following. In cal- 
culating the correlation function of the DC conductances 
at different energies, the energy difference AE enters for- 
mally as a finite frequency. This was first pointed out 
by Lee, Stone, and Fukuyama |^ in their diagrammatic 
evaluation of F in diffusive metals. The easiest way to 
see that ^ must enter the scaling function is to consider 
AE — 0, in which case, Eq. (0) gives the expected result 
@: F{Ep) = {5g^(Ep)) ^ r[L/£^{Ep)]. Moreover, both 
£,{Ep -\- AE) and £,{Ep) must enter as scaling arguments 
in because the range of AE wherein ^{Ep -\- AE) ~ 
^{Ep) is given by AE < const x ^{Ep)-'-^+^/''\ which is 
very small in the critical regime and vanishes much faster 
than the transition width W. 

Eq. (^) shows that in general, in the critical regime 
of the QHT, both Ep and AE enter the energy correla- 
tion function of the conductance. More important is the 
dual-role played by the Fermi energy difference. Writing 
£,{E) \E - E*]-" and ~ \AE\-^/^ with z the dy- 
namical scaling exponent, and setting one of the Fermi 
energy Ep = E* = 0, Eq. (0) becomes. 



F{AE) = T (\AE\L^/'', \AE\L' 



(3) 



Eq. (H) clearly shows that AE is a coupling constant 
conjugate to the static correlation (localization) length, 
and at the same time, a quantity analogous to a finite 
frequency conjugate to the length scale determined by 
the dynamical scaling exponent z. As a result, both the 
static and the dynamic critical properties enter the DC 
conductance correlation function. In general one expects 



that Harris criteria vz > \ holds. The two scaling ar- 
guments in Eq. compete and the correlation function 
must show a novel crossover from the regime dominated 
by static [AE ~ L~^/") fluctuations at large AE to that 
dominated by dynamic [AE ~ -^^^) fluctuations at small 
AE. Consequently, one expects the energy correlation 
range Ec to interpolate between Ec L^^^" at large 
AE and Ec ~ L"^ at smafl A^;. 

We next present a direct numerical calculation of the 
conductance correlation function F in Eq. (|^) for an inte- 
ger QHT in which the effects of electron-electron interac- 
tions are not considered |Tl| . In this theoretical noninter- 
acting analog of the true integer QHT in real materials, 
it is known that ~ 2.3 [|o| and the dynamical exponent 
z = 2. That z — 2 comes from the energy level spacing in 
a d = 2 noninteracting electron system and is consistent 
with the (anomalous) diffusive dynamics known at the 
noninteracting integer QHT [O. We will demonstrate 
that Ec indeed decays as L~^'^eLt large AE and as 
at small AE. Remarkably, the crossover region between 
the latter two behaviors is rather broad in AE over which 
we find Ec ^ L~^. 

For convenience, we choose to describe the transport 
in the integer quantum Hall regime using the Chalker- 
Coddington network model jig , !!?! . The latter is a 
square-lattice of potential saddle points (nodes) where 
quantum tunnelings between the edge states of the Hall 
droplets take place. With a choice of gauge the 
transfer matrix at each node is given by. 



/ cosh^ 
I sinh i 



sinh0 
cosh0 



(4) 



with a single real parameter 9. We have explicitly verified 
that introducing randomness in 6 does not change any 
of our results near the transition in agreement with the 
results of Ref. [|l^. Away from the nodes, the edge elec- 
trons move along the links (equipotential contours) with 
a fixed chirality set by the direction of the magnetic field 
and accumulate random Bohm-Aharonov phases. Note 
that besides the distribution of these random link phases, 
in Eq. (^ is the only parameter of the network. Chang- 
ing 9 amounts to varying the Fermi energy Ep across the 
QHT. In the rest of the paper, we will present results in 
terms of Ep which is related to the parameter 9 by the 
choice ol9 = 9c exp{Ep - E*) and sinh 9c ^ 1 |l7|. 

We have performed large scale numerical calculations 
of the two-terminal conductance. To this end, two semi- 
infinite ideal leads are attached to the left and right ends 
of the disordered network and periodic or open 

boundary conditions are applied in the transverse di- 
rection. Let us consider disordered networks having L 
columns of nodes and W/2 channels. Under such set- 
tings, the two-terminal conductance of a given sample 
with a fixed disorder realization is given by the Landauer 
formula 



giE) = -^Tr[ttt], 



(5) 
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where t is the {W/2) x {W/2) transmission matrix. In the 
transfer matrix approach, it is convenient to express g in 
terms of the (W x W) transfer matrix T since the 
latter is multiphcative across the L cohimns of scattering 
nodes in the network. Defining the ordered eigenvahies 
({Ai}) of the symplcctic matrix, H — T^T, by Xi — 
exp(27i) for i = 1, . . . , W, Eq. (||) can be written as. 



w/2 



1 



^ cosh' ME)) 



(6) 



Thus the calculation of the two-terminal conductance is 
transformed into that of the eigenvalues of the transfer 
matrix product H. It is known that constructing T^T 
by direct matrix multiplications is numerically unstable 
when the system size is large. We use here the stable 
numerical algorithm developed recently for large scale 
conductance calculations |^. The details of this algo- 
rithm have been discussed in Refs. |p|,|l5|. The basic idea 
is to maintain the stability of matrix multiplications us- 
ing the method of matrix UDR-decomposition, and to 
extract the eigenvalues using the method of orthonor- 
mal projection. Specifically, one can show following a 
sequence of UDR-decompositions, the n-th power of H , 
can be written as i?" — UnDnRn- In the limit of large n, 
typically less than 15, (i) t/„ is a unitary matrix of which 
the columns converge to the eigenvectors of H; (ii) Z3„ 
is a diagonal matrix and the eigenvalues of H is given 
by DnD^^^; and (iii) i?„ converges to a limiting right 
triangular matrix with unity on the diagonals. We next 
present the numerical results obtained using this algo- 
rithm on networks with L — W, and L up to 128 in units 
of the lattice spacing. 



calculated from Eq. for 49,000 disorder realizations 
at L = 128. It shows a remarkable skewed log-normal 
behavior as a result of the sharp fall off of P[gc] close to 
le^/h. 

In Fig. 2, we plot the conductance as a function of the 
Fermi energy in a typical sample for four different sample 
sizes. It is important to note that these reproducible 
spectra exhibit remarkably smooth oscillations with well 
defined peaks and valleys in the transition regime. More 
oscillations appear with increasing system size. Since the 
transition width shrinks as L~^/^ , the typical spacing 
between the peaks and the valleys, i.e. the correlation 
range Ec, must decrease systematically with increasing 
L faster than as L~'^/" . These features are very different 
from the conductance fluctuation spectrum in diffusive 
metals of mesoscopic dimensions. 




FIG. 2. Conductance fluctuation spectra in four typi- 
cal samples of L = 48, 64, 96, 128. 
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FIG. 1. Critical conductance distribution showing the 
skewed log-normal behavior (solid-line). 

The critical conductance fluctuations at the QHT have 
been studied recently 1^,1). The critical conductance 
Qc was found to be broadly distributed between and 
le'^/h with log- normal characteristics of the central mo- 
ments. Fig. 1 shows the distribution function of \og{gc) 



From an ensemble of up to 50,000 samples for each 
L, having individual fluctuation spectrum exemplified 
in Fig. 2, we calculated the conductance correlation 
function, F{AE), defined in Eq. (|). In Fig. 3, we 
plot the normalized conductance correlation function by 
the variance of the critical conductance, Fn(Ai?) = 
F{AE)/F{0), as a function of AE for different system 
sizes. The half widths Ai?i/2 of each correlation func- 
tion curve can be determined to extract the correlation 
range Ec versus L. Doing so, we obtained Ec ^ L~^. To 
understand this rather surprising result, we now perform 
a scaling analysis of the correlation function in Eq. (|^). 
Notice that the scaling function has two arguments orig- 
inated from two different length scales, ^ ~ |Ai?|^'' and 
|Ai?|^^/^. The scaling function should be domi- 
nated by the smaller one of the two. However, both of 
them are controlled by the same parameter, AE, i.e. the 
difference in the Fermi energy around the critical point. 
For the noninteracting QHT, v ~ 2.3 and 1/z = 0.5. 
Thus Lu; < C for AE < 1 whereas ^ < for AE > 1. 
We therefore introduce a new length scale X{AE), which 
interpolates between £^{AE) and L^{AE) as a function 
of AE, and rewrite Eq. (||) as 
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F{AE) = T{X{AE)/L). 



(7) 
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FIG. 3. The normalized conductance correlation func- 
tion versus the Fermi energy difference. The inset shows 
the scaling plot of the function in Eq. (Q). 

The complete behavior of X{AE) is obtained by de- 
manding that all finite size data of at different L and 
AE in Fig. 3 collapse onto a single scaling curve when 
F^ is plotted vs X/L. Such a scaling plot is shown in 
the inset of Fig. 3 where the curve represents the scal- 
ing function T in Eq. (|^). To our knowledge this is the 
first demonstration of the scaling behavior of conduc- 
tance correlations near a quantum phase transition. 
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FIG. 4. The crossover length X{AE) versus AE in- 
terpolating between static and dynamic fluctuation dom- 
inated asymptotics. 

The obtained X{AE) shown in Fig. 4. Indeed, A inter- 
polates between the correct asymptotic behaviors dom- 
inated by the static and the dynamical critical fluctua- 
tions: (i) For small AE, X ~ |Ai?|^^/^ and consequently, 
the correlation energy Ec ^ L^^. (ii) For large AE, 
X ~ lAEl'", implying Ec ^ i"^/''. (in) As Fig. 4 shows, 
the crossover regime between the asymptotic limits is 



very broad. Remarkably, over almost the entire crossover 
region, A exhibits a well defined power-law. 



A - |A£:| 



-i/C 



C = 0.96 ±0.05. 



(8) 



As a result, the correlation energy Ec ^ L^^ over this 
broad region. This is consistent with the conclusion 
drawn from the half- width analysis of Fig. 3, and pro- 
vides the microscopic mechanism by which such an un- 
usual phenomenon takes place. 

In summary, we have studied the energy correlation 
function of the conductance, a central quantity in meso- 
scopic physics, in the integer QHE. Although we have 
focused on the QHT in our numerical calculations, the 
physics discussed here is quite general and pertains to 
the critical regimes of quantum phase transitions that 
are driven by the location of the Fermi energy such as 
the metal-insulator transitions in disordered electronic 
systems. 
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